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Results from the Columbia lattice group study of the QCD finite temperature phase transition with dynamical 
domain wall fermions on 16'^ x 4 lattices are presented. These results include an investigation of the U{1) axial 
symmetry breaking above but close to the transition, the use of zero temperature calculations that set the scale 
at the transition and preliminary measurements close to the transition. 



1. Introduction 

Domain wall fermions (DWF) Q provide a way 
to separate the approach to the continuum limit 
from the approach to the chiral limit. At fixed 
lattice spacing the chiral limit is approached as 
the size Lg of the "fifth" direction is increased. 
Furthermore, since this regulator does not break 
flavor symmetry and since the computing cost is 
only linear in Ls it may be possible to simulate 
two flavor QCD with three light pions close to the 
flnite temperature transition using today's super- 
computers. 

Preliminary results of QCD thermodynamic 
studies with DWF were presented in ||^. Results 
from simulations on 8^ x 4 lattices indicated the 
presence of a low temperature phase with spon- 
taneous chiral symmetry breaking but intact fla- 
vor symmetry and a high temperature phase with 
the full SU{2) X SU{2) chiral flavor symmetry. 
Also, preliminary results on 16"^ x 4 lattices just 
above the deconflning transition indicated a dif- 
ference in the tt and S susceptibilities and screen- 
ing masses. This was preliminary evidence for 
the anomalous breaking of the C/a(1) symmetry 
above the transition. 

Here the above preliminary results are ex- 
tended in two directions: 1) The investigation of 
the f/yi(l) axial symmetry breaking is completed 
with results at two values of the coupling close 
to the transition. 2) A full study of the transi- 
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tion region on 16^ x 4 lattices is performed and is 
accompanied by scale setting calculations at the 
transition point to determine the transition tem- 
perature and the pion mass. Some of these results 
have already been reported in Q . Also see 

The DWF action used in this work is as in 
with the modiflcations as in 0. Some details on 
the numerical methods can be found in |^ . For a 
review on DWF see M and references therein. 



2. The ?7a(1) symmetry above Tc 

An important question concerns the role of 
the ?7^(1) anomaly above the finite temperature 
QCD transition. Earlier attempts to address this 
question ||Jl^ using staggered fermions could 
not produce conclusive results because staggered 
fermions break the Ua{^) symmetry already at 
the classical level. On the other hand, DWF 
bring this problem under control. In particular at 
Ls = oo massless DWF do not break the Ua{^) 
symmetry at the classical level. Also, the DWF 
Dirac operator can have exact zero modes [ pl| , 
the necessary property for anomalous breaking. 
In [l^ it was shown that this property is main- 
tained to a good approximation for fairly small 
masses and for Lg as small as w 10. 

The difference of the susceptibilities x of the tt 
and S is used as a measure of anomalous symme- 
try breaking. This difference was measured as a 
function of the bare quark mass m/onal6x4 lat- 
tice at P = 5.45 and (3 = 5.40 {pc « 5.325). The 
results are presented in figure 1. The lines are fits 
to co + C2m'j and have x^/d.o.f. « 1. The absence 
of a linear term indicates that for = 16 the chi- 



2 



ral symmetry is effectively restored. The ruf = 
extrapolated values are 0.26(6) at /3 = 5.40 and 
0.08(3) at /3 = 5.45. Although both are not zero 
by a statistically significant amount their value 
is small when compared with x-^ and xs which 
are ~ 8. Universality arguments require that if 
the QCD phase transition is to be second order, 
the anomalous U{1)a must be broken. It is an 
open question as to whether the small size of the 
U{1)a symmetry breaking seen here is sufficient 
to support this theoretical prediction that the 
two-flavor QCD phase transition is second order 
]l3| |. For quenched QCD thermodynamic stud- 
ies with DWF see , and with overlap 
fermions see p6| . 
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Figure 1. 16* x 4, = 16 and mg = 1.9. The 
fits are to cq -l-C2mj and the stars are the irif = 
extrapolated values. 

3. The transition 

The transition region (5.2 < /? < 5.4) is studied 
on 16^ X 4 lattices using to/ = 0.02, a domain wall 
height Too = 1.9 and Ls — 24. These values were 
chosen based on previous exploratory studies on 
smaller lattices [§,§1. In particular toq = 1.9 was 
found to produce two flavor physics and Lg = 24 
to be large enough to keep the residual chiral sym- 
metry breaking effects on the chiral condensate to 
the few percent level above the transition and to 
« 15% level below the transition. The results are 
shown in figure 2. The gauge part of the action 
is a standard Wilson plaquette action. 

In order to set the scale a simulation ||2^] on 
an 8^ X 32 lattice was done at /? = 5.325 corre- 



sponding to the middle of the transition region . 
It was found that in lattice units nip = 1.18(3) 
and TO^ = 0.654(3). This gives a value for 
the critical temperature Tc = 163(4) MeV and 
TOtt = 427(11) MeV. The critical temperature 
is in agreement with results obtained from other 
fermion regulators |^]. The pion mass is clearly 
too heavy to be able to extract useful information 
regarding the order of the transition. More so- 
phisticated studies |l^ indicate that the residual 
chiral symmetry breaking effects are much larger 
than expected and in order to obtain a physical 
pion mass Ls ~ 100 may be needed. 
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Figure 2. (V'V') vs. /3 with mo = 1-9, ruf = 0.02 
and Ls = 24. The circles correspond to 16"^ x 4 
lattice with an ordered initial configuration (ic), 
the crosses to a 16^ x 4 with disordered ic, the 
diamond to an 8'^ x 4 with ordered ic, and the 
cross to an S'^ x 4 with disordered ic. The gauge 
part of the action is a Wilson plaquette action. 

Large residual chiral symmetry breaking effects 
have been related to the "roughness" of the gauge 
field configurations ||ll|,^,^,^. Therefore a pos- 
sible "cure" could be to use improved gauge ac- 
tions. This idea was explored in zero tempera- 
ture quenched QCD by using an Iwasaki improved 
gauge action. Dramatic improvements were seen 
at couplings that correspond to A^t = 4 finite tem- 
perature quenched QCD just above the transition 
poll . These results prompted a study of the dy- 
namical QCD finite temperature transition using 
Iwasaki improved gauge action with ci = —0.331. 
The results are presented in figure 3. All param- 
eters are the same as in figure 2. 
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Figure 3. Same as in figure 2 but witli an 
Iwasaki improved gauge action with ci — —0.331. 

Again, in order to set the scale a simulation 
pO[ on an 8^ x 32 lattice was done at /3 = 1.9 cor- 
responding to the middle of the transition region. 
It was found that in lattice units nip — 1.163(21) 
and m^r = 0.604(3). This gives a value for 
the critical temperature Tc = 166(3) MeV and 
rriTr = 400(7) MeV. The critical temperature is 
in agreement with the Wilson plaquette action re- 
sults. However, the pion mass did not get lighter 
by a significant amount. 

4. Conclusions 

The difference of the susceptibilities of the tt 
and S was measured at two couplings just above 
the finite temperature transition. This difference 
was found to be non-zero by a statistically signif- 
icant amount indicating that Ua{^) remains bro- 
ken above but close to the transition. However, 
the value of the difference of the susceptibilities is 
much smaller than the susceptibilities themselves 
(by a factor of « 30 — 100 for each point) suggest- 
ing that the size of C^a(1) is indeed rather small at 
these couplings. It is an open question as to what 
effect this may have to the order of the transition. 

The transition region was studied by dynami- 
cal simulations on 16'^ x 4 lattices with Ls ~ 24 
and the critical coupling was located. At that 
coupling zero temperature simulations were done 
on 8^ X 32 lattices in order to set the scale. 
For the gauge part of the action two different 



types were used. A standard plaquette action 
gave critical temperature Tc — 163(4) MeV and 
m^r = 427(11) MeV. In an effort to reduce 
the pion mass an Iwasaki gauge action was also 
used. For that action = 166(3) MeV and 
m^r = 400(7) MeV. Unfortunately the pion mass 
is still too heavy to allow a study of the order of 
the transition. Larger values of Lg or alternative 
improvement techniques must be employed in or- 
der to reduce the pion mass to physical values 
and still be able to study the transition with the 
available computing resources. 

All calculations were done on the 400 Gflops 
QCDSP machine at Columbia University. 
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